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THE  OBSERVATIONS  OF  TURBULENT  BOUNDARY  LAYER  COHERENT  STRUCTURES  IN  AN  ADVERSE 
GRADIENT  FLOW 


Lian  Jixiang 

Beijing  Aeronautical  Engineering  College 

The  obterrationi  were  made  by  hydrogen  babble  method  is  a  water  chanael  6.S  m 
long  by  0.4  m  wide  and  0.3  m  deep.  A  flat  plate  6  m  long  was  set  along  one  side  of 
tbe  ehanneL  After  the  obaervations  of  TBL  atmetarea  is  aero  pressure  gradient  flow, 
a  cursed  wall  waa  aet  opposite  to  this  plate  to  fonn  a  two-dimenaioual  convergent- 
divergent  passage  as  shown  in  Fig.  1.  There  were  snuU  separated  bubbles  appearing 
intermittently  at  x  ^  200  to  300  mm  but  no  separation  in  the  downstream  part.  Plate  1 
is  a  plane  view  of  TBL  on  this  fist  piste  in  aero  pressure  gradient  flow  taken  at 
y*  ■«  8.4,  St  B  1200.  Plate  2  to  6  are  tbe  plane  views  and  tbe  side  views  of  TBL  on 
this  same  flat  plate  in  the  divergent  passage  at  x  »  100  to  300  mm  and  y  1  to  4  mm. 
In  this  adverse  pressure  gradient  flow  tbe  low  speed  spots  (regions  where  tbe  hy¬ 
drogen  bubble  lines  are  dense)  are  as  wide  as  tbe  high  speed  spots  as  shown  in  Piste 
2,  3,  4,  while  in  tbe  sero  pressure  gradient  flow  tbe  low  speed  regions  are  narrow  streaks 
(nanally  called  low  speed  streaks)  as  shown  in  Plate  L  In  tbe  interior  of  a  low  speed 
spot  (also  high  speed  spot)  the  hydrogen  bubble  lines  were  uniform  and  smooth,  there 
waa  no  trace  of  chaotic  flow.  Generally  tbe  chaotic  flow  appeared  only  at  one  of  tbe 
two  side  edges  of  each  low  speed  spot,  at  the  chaotic  side  tbe  hydrogen  bubble  lines 
concentrated  into  a  narrow  dense  white  long  streak  with  distorted  non-smooth  edges 
and  hydrogen  babbles  were  diffused  from  there,  while  at  the  opposite  side  the  hydrogen 
babble  lines  rmained  smooth-  Tbe  different  flow  features  of  the  two  side  edges  of 
each  low  speed  spot  were  ahown  clearly  in  Plate  2,  3,  A  These  dense  white  long 
•treaka  extended  far  down-stream  and  were  nmilar  to  tboee  streaks  in  sero  pressure 
gradient  Bows.  The  aide  views  show  many  large  transverse  vortices.  In  Plate  5  fluid 
at  *A*  was  raised  up  and  at  ‘C’  waa  pushed  down  to  wall  by  a  vortex.  In  Plate  6  tbe 
fluid  at  'A’  was  farther  from  tbe  babble  wire  than  the  free  stream  part  of  tbe  same 
bjrdrogea  babble  line,  therefore  it  abould  have  flowed  taster  than  the  free  stream. 


This  paper  was  received  in  April  1982. 


I .  Preface 


The  coherent  structure  of  turbulent  flow  is  an  active  branch  of  study 
originated  in  the  past  (more  than)  ten  years  -in  the  study  of  turbulent  flow;  the 
branch  plays  a  key  role  in  understanding  the  origination  regime  of  turbulent  flow. 
There  are  quite  a  few  papers  comprehensively  assessing  and  describing  the  coherent 
structure  [1-4]  of  the  boundary  layer.  There  have  been  many  papers  studying  the 
coherent  structure  of  the  turbulent  boundary  layer  with  zero  pressure  gradient; 
however,  there  are  few  papers  dealing  with  the  adverse  pressure  gradient.  Flow 
of  the  adverse  pressure  gradient  and  separation  of  the  boundary  layer  caused  by 
adverse  pressure  are  very  important  flow  phenomena,  related  (for  example)  to 
stall  of  aircraft  wings.  Kline  et  al.  [5]  and  Simpson  et  al .  [6]  conducted 
observations  of  the  adverse  pressure  zone  on  the  turbulent  boundary  layer; 
however,  there  have  been  no  reports  of  graphs  of  coherent  structure.  This 
experiment  was  conducted  on  a  two-dimensional  diffusion  sector  in  a  water 
channel  using  the  hydrogen  bubble  method.  The  paper  emphasizes  the  description 
of  observation  results  in  the  maximum  gradient  zone  of  adverse  pressure;  the 
boundary  layer  of  the  zone  is  close  to  separation.  At  the  beginning  of  the 
diffusion  sector  used  by  Simpson  et  al.  [6],  the  diffusion  angle  is  close  to 
zero.  The  diffusion  angle  along  the  flow  direction  gradually  increases;  the 
boundary  layer  is  entirely  separated  in  the  middle  portion  of  the  diffusion 
sector.  In  this  experiment,  the  diffusion  angle  is  the  greatest  at  the  begin¬ 
ning  of  the  diffusion  sector;  the  angle  gradually  decreases  along  the  flow 
direction  (refer  to  Fig.  1).  During  experiments,  the  diffusion  angle  should 
increase  as  much  as  possible,  so  that  the  boundary  layer  is  close  to  separation 
at  the  beginning  of  the  diffusion  sector.  Since  the  boundary  layer  at  the 
beginning  sector  is  thinnest,  a  relatively  high  gradient  of  adverse  pressure 
is  feasible;  thus,  no  separation  occurs  at  the  boundary  layer  at  the  end  of 
the  diffusion  sector.  Hence,  the  flow  field  is  relatively  stable,  providing 
convenience  in  observation.  The  observed  coherent  structure  is  significantly 
different  from  the  coherent  structure  of  a  turbulent  boundary  layer  without 
adverse  pressure  gradient. 


Fig.  1.  Flow  speed  of  free  stream  and 
boundary  layer  thickness  in  the  diffusion 
sector  and  the  shape  of  the  sector. 

II.  Experimental  Equipment  and  Method 

The  water  channel  is  6.8  meters  long  and  0.4  meter  wide  and  high;  the 
channel  is  made  of  organic  glass.  Near  the  lateral  wall  of  the  water  channel, 
an  organic-glass  plate  (6  m  long,  0.3  m  wide  and  25  mm  thick)  was  erected. 

First,  observations  on  turbulent  boundary  layer  without  pressure  gradient 
were  conducted  on  the  plate.  Then,  a  curved  wall  was  erected  near  another 
lateral  wall  of  the  channel.  A  two-dimensional  convergent -divergent  passage 
is  thus  formed  between  the  plate  and  the  curved  wall.  The  width  at  the  beginning 
of  the  divergent  sector  is  0.15  m;  at  a  distance  of  4.5  m  from  the  front 
edge  of  the  flat  plate,  the  diffusion  sector  is  1  m  long.  At  the  end,  the 
diffusion  sector  is  0.32  m  wide.  Refer  to  Fig.  1  for  the  shape  of  the  diffusion 
sector;  in  the  figure,  x  is  the  distance  from  the  cross  section  of  the  inlet 
of  the  diffusion  sector,  y  is  the  distance  from  the  surface  of  the  flate-plate 
wall . 

The  hydrogen  bubble  method  used  can  give  a  graph  of  turbulent  structure. 

A  platinum  wire,  25  in  diameter,  as  the  hydrogen-bubble-producing  filament, 

was  placed  in  water  for  observation.  The^ cathode  of  a  power  source  is  connected 
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with  the  platinum  filament,  and  the  anode  is  placed  in  water  in  a  non-inter¬ 
ference  flow  field;  the  power  source  is  square-wave  pulsed  electricity.  Its 
voltage  and  frequency  can  be  adjusted.  During  experimentation,  the  frequency 
is  fixed  at  9.1  Hz.  According  to  the  flow  speed  of  the  observation  zone  and 
the  length  of  the  platinum  filament,  the  voltage  can  be  determined;  in  this 
experiment,  the  voltage  used  was  15  to  150  V.  A  1  KW  projection  lamp  was  used 
for  illumination.  Under  the  action  of  pulsed  electricity  on  the  platinum 
filament,  hydrogen  bubbles  were  formed  in  streaks.  Each  pulse  produces  a  streak 
of  hydrogen  bubbles,  called  the  hydrogen  bubble  versus  time  line.  If  the 
platinum  filament  is  marked  with  coated  paint  or  a  pressure  trace  with  eqaal 
intervals,  the  equidistant  markings  also  form  hydrogen-bubble  streaks;  this  is 
called  the  hydrogen  bubble  versus  time  dyed  line.  All  these  methods  were  used 
in  the  experimentation. 

In  this  experiment,  the  water  depth  of  the  water  channel  was  0.3  m;  the 
flow  speed  at  the  inlet  of  the  diffusion  sector  was  19.2  cm/sec.  The  distribu¬ 
tion  of  flow  speeds  in  the  diffusion  sector  is  shown  in  Fig.  1.  The  observation 
proceeded  on  the  flat-plate  side  in  the  diffusion  sector;  the  front  edge  of  the 
flat  plate  is  the  turning  axis.  The  platinum  filament  was  made  to  be  perpendic¬ 
ular  to  the  flat  plate;  the  formed  hydrogen-bubble  streaks  were  generally 
situated  above  plane  xy.  The  observation  from  the  lateral  surface  (along  the 
X  axis)  is  called  the  lateral -surface  observation;  structures  like  a  vortex¬ 
shaped  movement  can  be  observed  within  the  boundary  layer.  The  platinum  filament 
was  made  parallel  to  the  flat  plate  and  perpendicular  to  the  free  stream;  the 
hydrogen-bubble  streaks  thus  formed  are  generally  parallel  to  the  wall  surface. 
The  observation  from  the  direction  perpendicular  to  the  flat  plate  is  called  the 
front-side  observation;  images  of  high-speed  spots,  low-speed  spots  and  streaks 
can  be  seen  at  the  bottom  layer  of  the  boundary  layer. 


In  the  side-view  photograph,  the  interval  of  two  hydrogen-bubble  streaks 
closest  to  the  platinum  filament  represents  the  instantaneous  flow  speed  at  the 
place.  By  using  a  microscope  coordinate  instrument,  the  distance  between  the  two 
front  edges  of  the  hydrogen-bubble  streaks  can  be  measured;  at  the  same  time, 
the  dimension  markings  on  the  measurement  negative  are  measured.  The  actual 
interval  between  two  hydrogen-bubble  streaks  is  then  obtained  by  proportional 


conversion;  the  actual  length  is  then  divided  by  the  time  interval  between  f  o 
pulses  to  obtain  the  instantaneous  flow  velocity  at  the  place.  From  tliis  method, 
the  pattern  of  instantaneous  velocity  of  the  boundary  layer  can  be  obtained;  then 
the  mean  velocity  pattern  can  be  derived.  For  details  of  this  method,  refer  to 
papers  [7  and  8].  By  using  this  method  in  experiments,  the  flow  velocity  U  of 
the  free  stream  in  the  diffusion  section,  and  the  distribution  of  thickness  5 
(of  the  boundary  layer)  along  the  x  axis  can  be  obtained,  in  addition,  a  hot- 
membrane  velocity-measuring  instrument  (of  50  diameter)  is  used  to  conduct 

the  same  measurement;  the  distribution  of  U  and  5  are  shown  in  Fig.  1. 

III.  E.xperimental  Results 

For  comparison,  first  observation  of  the  coherent  structure  of  a  turbulent 
boundary  layer  of  zero  pressure  gradient  is  conducted.  Then,  observation  is 
conducted  at  the  same  flat  plate  for  the  gradient  of  adverse  pressure  caused  by 
the  installation  of  a  curved  wall;  the  obtained  coherent  structural  graph  of  the 
turbulent  flow  is  significantly  different  from  that  of  the  zero  pressure  gradient. 
In  the  sector  from  x=100  mm  to  x=300  mm,  the  gradient  of  adverse  pressure  is 
greatest.  There  are  some  outstanding  characteristics  for  the  turbulent -flow 
structure.  In  the  following,  the  description  is  emphasized  of  the  observation 
results  of  this  sector  (simply  called  the  beginning  sector) . 

From  a  series  of  front-  and  side-view  graphs  observed  within  the  turbulent 
boundary  layer  of  a  flat  plate  of  zero  pressure  gradient,  the  results  are  general¬ 
ly  consistent  with  the  results  obtained  by  Kline  et  al.  [5]  and  Kim  et  al.  [7]. 
Photograph  1  of  plate  1  is  a  front-view  graph  of  the  experiment;  the  photograph 
shows  the  major  structural  images  of  fast-speed  spots  (within  the  zone  of 
sparse  hydrogen-bubble-time  streaks),  low-speed  spots  (within  the  zone  of  dense 
hydrogen-bubble-time  streaks),  and  streaks  (long  white  streaks  formed  downstream 
of  the  block  of  spots  by  hydrogen  bubbles,  generally  along  the  direction  of  the 
free  stream).  In  photograph  1  of  plate  I,  hydrogen-bubble-time  streaks  are  used; 
these  streaks  are  different  from  the  hydrogen-bubble-time  dyed  lines  (equal - 
interval  coated  paint  markings  on  the  platinum  filament)  used  by  Kline  et  al.  [5]. 
In  appearance,  these  two  graphs  are  different.  However,  by  examining  them 

closely,  the  characteristics  of  the  two,  in  that  they  show  the  fast-speed  spots, 
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at  approximately  100),  are  consistent. 


However,  in  the  paper  the  images  of 


fast-speed,  low-speed  spots  and  streaks  shown  by  hydrogen-bubble-time  streaks  are 


clearer.  The  streak  graphs  shown  by  the  two-cases  are  generally  consistent. 


The  x-direction  width  of  low-speed  spots  in  a  zero-pressure  gradient  is  very 
narrow.  Even  at  places  very  close  to  the  wall  surface  (y^=3.5  as  achieved  by  the 
author  and  his  colleagues),  the  width  of  low-speed  spots  is  much  smaller  than 
that  of  fast-speed  spots.  Therefore,  the  shape  of  low-speed  spots  is  generally 
narrow  streaks  along  the  stream  direction.  At  the  beginning  sector  of  the 
experiment,  the  x-direction  width  of  low-speed  spots  is  close  to  or  exceeds 
the  fast-spot  width,  as  shown  by  photographs  2  through  4  in  plate  I.  Therefore, 
all  graphs  of  low- speed  and  high-speed  spots  are  spots  of  large  area.  Compared 
to  the  zero  pressure  gradient,  the  hydrogen- bubble  streaks  within  low-speed 
spots  are  denser,  and  the  difference  in  densities  of  hydrogen-bubble  streaks  of 
the  low-speed  spots  and  fast-speed  spots  is  greater.  This  is  because  the  flow 
speed  within  low- speed  spots  is  lower.  The  low-speed  spots  become  wider  and 
thicker,  and  the  flow  velocity  within  the  spots  is  lower;  these  are  characteris¬ 
tics  for  a  boundary  layer  close  to  separation.  This  further  lowers  the  flow  speed 
within  the  low- speed  spots;  even  the  flow  speeds  become  negative  values,  thus 
forming  reverse-flow  spots  (a  divergent  bubble) .  In  the  middle  portion  of 
photograph  4  of  plate  I,  there  are  a  few  hydrogen-bubbles  flowing  to  the  upper 
stream  of  the  platinum  filament;  this  is  a  reverse-flow  spot.  In  the  beginning 
sector  within  the  zone  of  x=200  mm  to  300  mm,  small  reverse-flow  spots  generally 
appear  intermittently,  revealing  the  closeness  of  divergence  of  the  boundary 
layer  at  the  place.  Therefore,  in  the  paper,  the  dimensionless  wall  surface  units 
(such  as  y*  and  z^)  are  not  used  to  indicate  the  experimental  data  in  the  adverse 
pressure  zone. 


Since  the  low-speed  spots  are  wider,  we  can  clearly  discriminate  the  flow 
images  inside  the  low-speed  spots  and  their  edge.  We  can  see  from  photographs  2 
through  4  in  plate  I  that  the  distribution  of  hydrogen-bubble  streaks  in  low-speed 
spots  is  relatively  even,  like  that  inside  the  high-speed  spots.  Intervals 
between  streaks  are  changed  moderately  along  the  Z  direction;  this  indicates  that 
the  flow  velocities  within  spots  are  relatively  even.  The  relatively  high 


velocity-gradient  appears  in  the  transitional  zone  at  the  junction  point  of 
low-speed  and  high-speed  spots.  From  photographs  2  through  4  in  plate  I, 
and  thousands  of  pictures  and  films  photographed  by  the  author  and  his  colleagues, 
we  can  see  that  for  two  lateral  edges  of  a  low-speed  spot,  generally  only  one 
side  is  rapidly  transferred  into  a  fast-speed  spot;  on  the  other  side,  it  is 
moderately  transferred  to  a  fa'  eed  spot.  At  the  lateral  edge  of 
moderate  transition,  the  inter  s  of  hydrogen-bubble  streaks  generally  change 
from  dense  to  sparse,  forming  dder  transition  belt.  The  curve  of  each 
hydrogen-bubble  streak  is  smoo  n  .th  its  edge  clearly  visible,  without  a  trace 
of  chaotic  flow.  At  the  side  of  fast  transition,  the  intervals  of  hydrogen- 
bubble  streaks  are  rapidly  becoming  smaller.  Not  far  from  downstream  of  the 
platinum  filament,  all  hydrogen- bubble  streaks  converge  densely  together,  forming 
a  long,  slender  and  dense  white  streak;  this  narrow  streak  separates  the  fast- 
speed  and  low-speed  spots;  the  transcendent  streak  has  a  very  high  velocity 
gradient.  This  white  streak  generally  extends  along  the  flow  direction  toward 
downstream,  as  indicated  by  "  "  on  photograph  2  of  plate  I;  the  shape  is  not 
much  different  from  the  streak  of  the  bottom  layer  of  a  zero  pressure  gradient. 

We  can  see  that  the  edge  of  these  streaks  is  not  smooth,  with  hydrogen  bubbles 
diverging  from  the  edge;  the  streaks  are  twisted  or  divergent  downstream.  These 
graphs  indicate  turbulent  motion.  We  can  also  see  by  comparing  these  graphs 
the  following:  the  slower  the  low-speed  spots  (in  other  words,  the  denser  the 
hydrogen-bubble  streaks),  the  more  chaotic  are  the  streaks  at  the  side.  By 
visual  observation,  direction  vortices  frequently  appear  at  the  side  edge  of 
the  rapid  transition  zone.  From  the  right  side  of  the  reverse-flow  spots  of 
photograph  4  of  plate  I,  the  streaks  are  spiral  shaped;  downstream,  there  is 
a  large  conic  vortex.  There  is  no  chaotic  flow  at  the  side  of  moderate 


By  observing  the  lateral  surface  of  the  beginning  sector,  in  the  majority 
of  photographs  there  are  transverse  vortex  motions,  but  very  few  vortices  along 
the  flow  direction.  This  phenomenon  is  different  from  that  of  a  cero  pressure 
gradient.  In  the  observation  results  by  Kim -et  al.  [7]  on  the  lateral  surface 
of  a  turbulent  boundary  layer  of  zero  pressure  gradient,  there  are  mostly  flow- 
direction  vortices,  estimated  as  two  thirds  of  all  vortices  and  one  third  of  the 
lateral -direct ion  vortices.  At  the  beginning  sector,  there  are  many  large 
lateral -direction  vortices;  it  can  be  seen  that  these  vortices  exert  an  influence 
on  the  flow  of  most  zones  of  the  local  boundary  layer.  In  the  second  hydrogen- 
bubble  streak  (the  hydrogen-bubble  streak  closest  to  the  platinum  filament  is  the 
first  one,  then  in  sequence  toward  the  downstream  in  the  order  second,  third, 

....)  of  photographs  2  through  4  of  plate  I,  the  vortex  center  is  about  5  mm 
from  the  wall  surface.  By  its  influence,  the  left  lower  terminal  of  the  second 
hydrogen- bubble  str'ak  rises  to  A;  its  right  side  is  lowered  to  C.  In  the 
photograph,  B  is  the  left  lower  terminal  of  the  third  hydrogen-bubble  streak. 
Without  the  influence  of  this  vortex,  it  should  be  close  to  the  wall  surface,  but 
actually  it  rises,  with  only  a  narrow  seam  remaining  with  the  sector  AC  of  the 
second  hydrogen-bubble  streak.  The  induction  of  a  lateral -direction  vortex 
reduces  the  flow  velocity  while  raising  the  flow  at  its  lower  rear  side.  In 
photograph  5,  there  is  a  gap  in  the  lower  terminal  of  the  first  hydrogen-bubble 
streak,  indicating  that  the  flow  velocity  is  reduced,  as  it  is  influenced  by  the 
preceding  vortex.  The  induction  of  a  lateral-direction  vortex  on  its  upper  front 
fluid  shows  relatively  clearly  in  photograph  6  of  plate  I.  In  the  photograph, 

A  is  situated  at  the  right  upper  side  of  a  lateral -direction  vortex.  At  the  place 
influenced  by  the  vortex  on  the  dyed  streaks  (the  black  strr iks  in  the  hydrogen- 
bubble  streaks  formed  by  pressure  traces  on  the  platir.jm  filament)  ,  the  hydrogen 
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raising  the  fluid.  This  function  of  the  lateral-direction  vorte.x  creates 
an  exchange  of  momentum  on  a  large  scale. 


IV.  Conclusions 


This  paper  presents  the  observation  results  of  a  two-dimensional  divergent 
sector  (the  sector  from  x=100  mm  to  x=300  mm)  in  a  water  channel.  This  sector 
is  the  largest  portion  of  the  adverse  pressure  gradient  of  the  pressure 
spreading  sector.  During  e.xperiments ,  a  small  number  of  divergent  bubbles 
intermittently  appear,  indicating  the  closeness  of  divergence.  However,  down¬ 
stream  from  this  sector,  there  are  no  divergent  bubbles.  Compared  to  the  zero 
pressure  gradient,  the  main  characteristics  of  the  structure  of  the  turbulent 
boundary  layer  are  wide  low-speed  spots  and  further  lower  flow  velocity,  more 
difference  from  the  flow  speed  of  high-speed  spots;  there  are  many  large 
transverse  vortices.  Since  the  low-speed  spots  are  wider,  the  flow  can  be 
clearly  discriminated  between  the  inside  and  the  edge  of  the  low-speed  spots. 
Within  the  low-speed  and  high-speed  spots,  the  flows  are  relatively  homogeneous 
and  steady;  at  the  same  time,  the  other  side  is  chaotic.  The  hydrogen-bubble 
streaks  of  the  steady  side  are  smooth  and  clear  with  a  wider  transitional  zone 
between  the  high-speed  spots.  In  the  chaotic  side,  the  hydrogen-bubble  streaks 
are  dense,  becoming  a  beam  of  a  dense  white  streak.  The  edge  of  the  streak  is 
not  smooth,  with  divergent  hydrogen  bubbles,  or  with  traces  of  the  chaotic 
motion,  such  as  twisted  shape  or  divergence.  The  streak  extends  far  downstream 
its  shape  is  similar  to  the  streak  of  a  zero  pressure  gradient.  In  the  quiet 
side,  generally  no  streaks  are  formed;  the  transverse  vortices  are  larger.  The 
flow  velocity  in  the  vicinity  of  some  vortices  is  apparently  greater  than  the 
flow  velocity  outside  the  local  boundary  layer. 
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THREE-DINENSIONAL  TURBOMACHINE  TRANSONIC  FLOW  EQUATIONS  AND  CHARACTERISTIC 
THEORY  EXPRESSED  WITH  RESPECT  TO  NON -ORTHOGONAL  CURVILINEAR  COORDINATES 


Lu  Wenqiang 

Institute  of  Engineering  Thermophysics,  Chinese  Academy  of  Sciences 


By  the  nse  of  tensor  ealeolna,  three-dimensional  torbomschine  transonic  flow  equa¬ 
tions  in  weah-conser-Tstive  form,  ^e-dependent  iterstive  equations  of  two  kinds  of 
fluid  surface  and  eharaeteristie  theory  expressed  with  respect  to  non-orthogonal 
curvilinear  coordinates  have  been  derived  in  this  paper. 


This  paper  uses  tensor  calculus  to  derive  the  three-dimensional  weak 
conservation  form  equations,  time  dependent  flow  surface  iterative  equations, 
and  characteristic  theory  of  non-orthogonal  curvilinear  coordinates.  These 
equations  and  theory  are  valuable  in  solving  the  three-dimensional  transonic 
flow  field  (with  shock  wave)  of  turbomachines,  in  addition  to  numerical 
analysis  of  circumferential-  and  diametral -direct ion  anomalies,  and  for  calcu¬ 
lating  the  trembling  vibration  unsteady  aerodynamic  force. 

1,  Three-dimensional  weak  conservation  form  equations  of  non-orthogonal 
curvilinear  coordinates 

The  non-orthogonal  curvilinear  coordinate  three-dimensional  flow  equations 
[2]  of  turbomachines  are 

This  paper  was  received  on  4  May  1982. 
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By  using  the  relationship  equation: 


-  V/  +  TV5  -  -  -1 V,  -  V  ilCHO*  -  («0’1; 

<1^,  -  Vi"w>oo.«.r. 

vri 

(2)  X  pV^  +  (1)  X  «■,;  D)  X  ^>/T  +  (1)  X 

O)  X  pTT  +  (1)  X  «/,;  (5)X^+(1)X/ 


We  then  derive  the  three-dimensional  weak  conservation  form  equations 
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2.  Time  dependent  flow  surface  iterative  equations 

The  author  and  his  colleagues  derived  fundamental  equations  (along  the  flow 

1  T  3  4 

surface)  on  the  curvilinear  coordinate  system  (x  ,  x* ,  x  ,  x  )  fixed  at  the 

rotating  shaft  of  turbomachines.  In  the  fundamental  equations,  x*^  =  t  and  the 

basic  vector  J.  =  ^/5x^. 

1 

3  3  1 

In  the  space,  there  is  an  unsteady  flow  surface  x  =x  (x  ,  x  ,  t) . 

For  the  normal  and  tangential  direction  dot  product  of  the  flow  surface: 

II  •  rfr  —  n,Jx'  —  -^=  ir'  —  0 

Tti  AnitA^^nn  *  ^ 


In  the  equation, 


f.i  “  Ji  ’  Jr,  inf'  ”  «» ;  —  It  2. 3, 1 
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(1)  Continuity  equations 


At  an  instant,  by  taking  any  element  of  the  flow  surface  with  the 
measurement  —  y/~*  dx'dx*;  and  by  prescribing  the  area  density  j  of 

the  flow  sector,  then  the  flow  mass  on  an  element  of  flow-sector  surface  is 

ads  m 

According  to  the  law  of  conservation  of  mass  and  Reynolds'  transmission 
theory,  we  obtain 


and  then 


~  +  Av^iaW)  -  0 


Since  the  aforementioned  computations  proceeded  at  a  flow  surface,  in  order 
to  discriminate  the  conventional  computation,  the  following  form  of  the  equation 


IS  written: 


3  (pr>/V )  S  (pr«'*V  *  ) 


In  the  equation,  the  thickness  (of  the  flow  sector)  is  r  —  p/a  —  •  )</x’ 

(2)  Momentum  equations 

In  order  to  cancel  out  the  partial  derivatives  in  the  x"’  direction,  the 
relationship  W-n=0  is  applied  in  addition  to  the  following  calculations: 

Ja - _ 

3x*  Ox*  Ox*  Ox*  Ox*  ■*  Or* 

Canceling  out  derivatives  along  the  direction  in  equations  (2),  (3)  and  (4), 
we  then  obtain: 


2y7(«^«*  - 

8t  1 3x‘  Ox*  J  Ox* 


^  +  -  -  uW) 

Bt  I  Ox*  Ox*  J  Ox* 
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(3)  Time  dependent  flow  surface  iterative  weak  conservation  form 

equations 


Processing  Eqs.  (71  and  (8)  as  in  section  1,  we  then  obtain  the  time 
dependent  flow  surface  iterative  weak  conservation  form  equations 
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In  the  equations,  the  steady  state  energy  equation  (dl/dt=0)  is  used  for  the 
energy  equation,  as  expressed  by  numerical  experiments,  thus  saving  time  in 
reaching  the  steady  state. 

3.  Characteristic  theory 


The  characteristic  theory  describes  the  objective  rule  of  message  exchange 
between  the  boundary  and  the  internal  flow  field.  By  using  the  boundary  condi¬ 
tions  established  by  the  characteristic  theory,  the  time  progressing  process 
of  numerical  solution  can  be  ensured  with  the  reliable  mathematical  basis  of 
physics. 


Equations  (7),  (8)  and  equations  of  unsteady,  cohesionless  and  adiabatic 
energy  equations 

Bt  B*‘  I  8*' i 


are  condensed  and  combined  into 


2  a5  •  gndmg  AT,  ^“1,2,  3,  4,  5 

Kxl 

The  aforementioned  equations  are  combined  linearly 


“  K»1  **•* 

If  we  assume  the  normal  direction  (of  the  characteristic  cone) 

A=A.e^  1»  2 

J 

Then  we  have  the  characteristic  equation 

(10) 

The  characteristic  compatibility  relationship  is 


(11) 


In  the  equation,  a  is  a  root  of  the  equation 


-  0 

(1)  Characteristic  equation 


(1-0 
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We  then  obtain  the  following  characteristic  types: 

2  2 

a.  0=^-W=0  streamline;  b.  D  =C  sonic-wave  envelope, 

(2)  The  characteristic  compatibility  relationship  of  the  non-orthogonal 
curvilinear  coordinate  system 

a.  If  we  substitute  D=0  into  Eq.  (12),  we  then  obtain  the  following 
flow  characteristic  compatibility  relationship: 

2  3 

The  solution  of  the  first  group:  ®  ^24^^^  '^34”®  momentum  equation 
combination; 

The  solution  of  the  second  group:  a^=0  energy  equation 

2  2 

b.  If  we  substitute  D  =C  into  Eq.  (12),  we  then  obtain 

«*  —  C’;  «*  —  —  pt-ZT Da»;  «•  —  —  pr-v/T Dl*; 

a*—0;  *>  — 

If  we  substitute  the  aforementioned  equations  into  Eq.  (11),  we  then 
obtain  the  wave  characteristic  compatibility  relationship; 
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By  using  the  characteristic  theory  at  the  boundary,  and  conducting  a 
similar  discussion  like  paper  [2],  we  obtain  the  corresponding  relationship. 
For  three-dimensional  equations  (6),  a  similar  method  can  be  used  to  derive 
the  characteristic  theory  of  non-orthogonal  curvilinear  coordinates,  and  to 
conduct  the  corresponding  discussion  at  boundaries. 


4.  Conclusion 


(1)  The  introduction  of  non-orthogonal  curvilinear  coordinates  can 
allow  us  to  select  the  most  suitable  and  practical  curved  lattice  for  a  comple.x 


shaped  blade  grid  and  different  types  of  flow  problems  in  order  to  raise  the 
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calculation  accuracy  and  convergent  speed.  This  compensates  shortcomings  in 
the  old  time-dependent  method  where  arbitrary  curve  lattices  cannot  be 
formed.  This  is  the  reason  that  in  the  early  1980s,  the  body-fitted  coordinate 
system  was  applied  in  several  transonic  numerical  methods  [5-5]  abroad.  .At 
the  appointed  revolving  surface,  we  conducted  preliminary  numerical  e.xperi- 
ments  [2].  Our  work  shows  that  the  application  of  non-orthogonal  curvilinear 
coordinates  in  the  time  dependent  method  is  successful. 

(2)  Ways  of  solving  the  three.-diraensional  transonic  flow  field  (with 
shock  wave)  of  turbomachines  had  been  preliminarily  studied  abroad.  For  example, 
in  [6]  which  applied  the  weak  conservation  form  equation  of  a  three-dimensional 
cylindrical  coordinate  system  by  selecting  the  MacCormack  time-splitting  format, 
in  the  space,  first  proceed  with  time  increments  along  the  r  direction,  and 
then  the  <P  direction,  and  finally  time  increments  along  the  z  direction.  The 
three-dimensioanl  transonic  flow  field  was  obtained  by  using  large  amounts  of 
computer  time.  In  manuscripts  [7]  and  [8],  the  time  dependent  iteration  is 
similar  to  Sj^-S2  (two  types  of  flow  surface)  of  the  cylindrical  coordinate 
system.  By  using  less  computer  time,  the  three-dimensional  approximate  flow 
field  of  the  cylindrical  coordinate  system  can  be  obtained  which  is  comparable  to 
experiments.  We  consider  that  the  strict  time  dependent  two-category  flow 
surface  iterative  method  can  have  time  increments  of  the  region  to  be  solved, 
approximating  the  restraint  of  the  flow  surface.  Obtaining  the  complete  three- 
dimensional  solution  of  the  entire  flow  field  is  completed  with  the  establishment 
of  the  steady-state  solution  of  various  groups  of  flow  surfaces,  and  the  iterative 
convergent  solution  of  two  groups  of  flow  surfaces.  This  allows  the  skillful 
application  of  time  increments  of  a  single  group  of  flow  surfaces,  and  an 
alternate  iteration  technique  of  two  groups  of  flow  surfaces.  Of  course,  the 
advantages  and  disadvantages  of  the  two  aforementioned  methods  of  solving  the 
three-dimensional  flow  field  are  subject  to  comparative  numerical  experiments. 

The  author  expresses  his  thanks  to  Professor  Wu  Zhonghua  for  his  counsel 
and  support  in  this  research  work. 
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of  photographs  2  through  4  of  plate  I,  the  vortex  center  is  about  5  mm 

from  the  wall  surface.  By  its  influence,  the  left  lower  terminal  of  the  second 

hydrogen-bubble  str  ak  rises  to  A;  its  right  side  is  lowered  to  C.  In  the 

photograph,  B  is  the  left  lower  terminal  of  the  third  hydrogen-bubble  streak. 

Without  the  influence  of  this  vortex,  it  should  be  close  to  tlie  wall  surface,  but 

actually  it  rises,  with  only  a  narrcw  seam  remaining  with  the  sector  AC  of  the 

second  hydrogen-bubble  streak.  The  induction  of  a  lateral -direct  ion  vortex 

reduces  the  flow  velocity  while  raising  the  flow  at  its  lower  rear  side.  In 

photograph  5,  there  is  a  gap  in  the  lower  terminal  of  the  first  hydrogen-bubble 

streak,  indicating  that  the  flow  velocity  is  reduced,  as  it  is  influenced  by  the 

preceding  vortex.  The  induction  of  a  lateral -direct ion  vortex  on  its  upper  front 

fluid  shows  relatively  clearly  in  photograph  6  of  plate  I.  In  the  photograph, 

A  is  situated  at  the  right  upper  side  of  a  lateral-direction  vortex.  At  the  place 

influenced  by  the  vortex  on  the  dyed  streaks  (the  black  strr iks  in  the  hydrogen- 

bubble  streaks  formed  by  pressure  traces  on  the  platinum  filament) ,  the  hydrogen 

bubble  streaks  have  turned  about  45°  along  the  clockwise  direction.  Ai  the 

place,  the  hydrogen-bubble  streaks  are  farther  from  the  platinum  filament 

than  the  portion  of  the  same  hydrogen-bubble  streak  outside  the  boundar)  layer. 

This  indicates  that  within  the  duration  (about  two  pulse  periods)  from  the 

platinum  filament,  the  flow  velocity  of  the  hydrogen-bubble  streak  near  the 

wall  is  greater  than  the  flow  velocity  of  the  free  stream.  Apparently,  this  is 

the  result  of  induction  by  the  vortex.  In  many  photographs,  the  lateral - 

direction  vortices  speed  up  the  fluid  at  its  upper  front  side  to  let  the  fluid 

descend,  and  reduce  the  velocity  of  the  fluid  in  the  rear  downward  side  by 
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